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Abstract 

We study extreme-ultraviolet emission line spectra derived from three-dimensional 
magnetohydrodynamic models of structures in the corona. In order to investigate 
the effects of increased magnetic activity at photospheric levels in a numerical ex- 
periment, a much higher magnetic flux density is applied at the photosphere as 
compared to the Sun. Thus, we can expect our results to highlight the differences 
between the Sun and more active, but still solar-like stars. We discuss signatures 
seen in extreme-ultraviolet emission lines synthesized from these models and com- 
pare them to observed signatures in the spatial distribution and temporal evolution 
of Doppler shifts in lines formed in the transition region and corona. This is of 
major interest to test the quality of the underlying magnetohydrodynamic model 
to heat the corona, i.e. currents in the corona driven by photospheric motions (flux 
braiding). 
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1 Introduction 



The mechanism of field line braiding as a heating mechanism of the so- 
lar corona has been discussed for a long time. However, three-dimensional 
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numerical models of the solar corona and solar-like stars based on magneto- 
hydrodynamics only became available in recent years (Gudiksen & Nordlund 



20021 . l2005al . l2005bl ). A small active regio n was described in a numerical ex- 



periment by Gudiksen & Nordlund (120021 ) resulting in a corona consisting of 
small loop-like structures. In order to investigate the general appearance of 
the synthesized cor ona, synthe tic spectra can be produced from the simulated 
data (Peter et al. 12004 . 120061 ) and can then be compared to spectroscopic 
observations of the Sun. 



In this work, we model the atmosphere above an active region in a three- 
dimensional box that extends from the photosphere to the corona. The lower 
boundary magnetic field consists of two parts, i.e. an active region magnetic 
field and a quiet Sun network magnetic field. The active region magnetic field 
is the same as in Gudiksen & Nordlund (120021 ). i.e. a spatially downscaled mag- 
netogram of Active Region 9114 observed by the Michelson Doppler Imager 
(MDI; Scherrer et al. Il995l ) onboard the Solar and Heliospheric Observatory 
(SOHO) on August 8, 2000. The downscaling was performed in order to fit 
the active region into the computational domain (50x50 Mm 2 horizontally). 
As the downscaling suppresses the small-scale network fields, a typical quiet 
Sun MDI high-resolution magnetogram was added. In order to highlight the 
interaction of the active region with the neighboring quiet Sun, the quiet Sun 
flux density was enhanced by a factor of four. By this a generic magnetogram 
of an active region is produced that allows us to study the interaction of 
magnetic structures on small scales. As a lower boundary condition a time- 
dependent v elocity field with properties close to so lar granulation (Gudiksen 
& Nordlund [2002]) is imposed. Peter et al. (120041 ) reproduced the observed 
variation of line shifts with temperature on the Sun from three-dimensional 
numerical models, especially the temperature of the maximum redshift. This 
temperature (logT ss5.3) compares well to the v alues found fo r activ e region 
and quiet Sun observations (Peter & Judge Il999l . Teriaca et al. Il999l ). There- 
fore, our numerical experiment can give a first indication of what to expect in 
a situation with a higher magnetic complexity. 

Active stars with convective envelopes show a complicated surface magnetic 
field distribution. High-resolution spectra are needed to separate the measured 
magnetic flux Bf, i nto the magnetic field density B and the magnetic filling 
factor /. Saar (119961 ) investigated the behavior of stellar surface magnetic fields 
of main sequence dwarf stars and found a correlation between the measured 
magnetic flux Bf, and the rotation period of a star. The stellar magnetic flux 
Bf, as well as the filling factor / of stars that rotate slower than the satu- 
ration threshold grows with larger rotation rate. The fluxes reach values up 
to 100 times larger than those observed for the Sun and i mply an increased 
heating of the upper layers of the stars' atmospheres. Saar (119961 ) claims that 
saturation occurs in the filling factor /, but not in the magnetic flux Bf. From 
this result one could conclude that a star rotating at the saturation threshold 
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is completely covered with magnetic field (/=1) and exhibits a much higher 
magnetic complexity compared to the Sun, where the magnetic flux is concen- 
trated in relatively small regions of strong fields. Doppler imaging and Zeeman 
Doppler imaging give evidence for s tarspots and stellar surface magnetic flux 
at high latitudes (Donati et al. Il999l). but also for spots and magnetic fields at 
low latitudes (e.g. Barnes et al. Il998l ) suggesting more complex surface mag- 
netic field topologies than found on the Sun. The setup of the model we use 
in this study to investigate properties of the synthesized emission lines can be 
considered as a step towards higher magnetic complexity and higher magnetic 
flux at surface levels, and thus as a step towards investigating the coronae of 
stars more active than our Sun. 



The temperatures and densities in the simulation box as well as the veloc- 
ity components and the magnetic field components are accessible at each grid 
point for every time-step of the simulation. The model temperatures and den- 
sities cover several orders of magnitude. The well-known rise in temperature 
and the drop in density from the photosphere to the corona are reproduced 
well in the box. We synthesize emission-line spectra from such coronal mod- 
els and investigate average properties of these spectra, such as intensities and 
Doppler shifts. 



2 MHD model and spectral analysis 



In this model the heating of the corona is based on the braiding of mag- 



netic fi eld lines due to motions in the photosphere, as first suggested by [Parker 



(119721 ). The braided magnetic field produces currents that are dissipa ted b y 



Ohmic heating. We are using the Pencil Code (Brandenburg & Dobler 120021 ) . 
a high-order finite difference code (sixth order in space and third order in 
time) for solving the MHD equations in a compressible medium. The three- 
dimensional magnetohydrodynamic model accounts properly for the energy 
balance, especially for heat conduction and radiative losses, allowing us to 
reliably synthesize the profiles of extreme-ultraviolet emission lines observ- 
able with current EUV spectrometers, i.e. the Solar Ultraviolet Measurement 
of Em itted Radiation (SUMER) instrument onboard SOHO (Wilhelm et al. 



19951 ) and the Extr eme- u ltraviolet Imaging Spectrometer (EIS) onboard Hin- 



ode (Culhane et al. 120071 ) . 



The computational box extends 50 Mm x 50 Mm in the horizontal direction 
and 30 Mm in the vertical direction covering the atmosphere from the photo- 
sphere to the corona with a 256 3 grid that is non-equidistant in the vertical 
direction. The vertical grid spacing goes down to 80 km in the low transition 
region (logT pa 4.3) around the maximum temperature gradient with a spac- 
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ing in temperature of 0.18 dex on a log T scale. For the calculation of the 
emission lines, the ve rtical spacing of t he grid is not sufficient. Therefore, a 



spatial interpolation (IPeter et all 120061 ) has to be applied. 



We determine the emissivities, intensities and Doppler velocities for a num- 
ber of extreme-ultraviolet emission lines (Table 1) from the temperatures and 
densities at each point of the interpolated grid in the box. Also shown in Table 
1 are the line formation temperatures of the respective line, i.e. the tempera- 
ture where the maximum of the contribution function is found. This is roughly 
comparable, but not identical to the temperature of the maximum ionization 
fraction. The lines have been selected for two reasons. First, they cover the 
range of temperatures in the transition region, i.e. from T fh 10 4 — 10 6 K. 
Second, the lines are observable with spectrometers such as SUMER/SOHO 
and EIS/Hinode, and are major lines in the bandpasses of the future Atmo- 
spheric Imaging Assembly (AIA) for the Solar Dynamics Observatory (SDO). 
The lines are optically thin, and predominantly excited by electron collisions. 
The emissivity of the transition from an upper level 2 to a lower level 1 is 
given by e = /iz/n 2 v4 2 i, where hv is the energy, n 2 is the number density of the 
upper level, and A 2 \ is the Einstein coefficient. To calculate the line emissivi- 



ties from the MHD model we employ the CHIANTI package (Dere et al. 11997 



Landi et al. 120061 ). For the calculation of the number density of the u pper 



level ionization equilibrium is assumed and the Mazzotta et al. (119981 ) ion- 
ization fractions are used. The validity of ionization equilibrium for modeling 
synthetic spectra from three-dimensi onal M HD models is a critical point and 



has been investigated by Peter et al. (120061 ). Their discussion of ionization and 



recombination times in comparison to the dynamic timescale, i.e. the time to 
cross the line formation temperature, showed that the ionization and recom- 
bination times in the coronal and transition region plasma are mostly smaller 
than the typical hydrodynamic timescales, indicating that the approximation 
of ionization equilibrium is not too bad. However, for future models the inclu- 
sion of non-equilbrium ionization is very important. Since we do not consider 
line ratios here, the results do not depend on the specific abundances used in 
this st udy. W e refer to the photospheric elemental abundances (Grevesse & 



Sauval ll998l ) in CHIANTI. The spectral profiles are assumed to be Gaussian 
with a line width corresponding to the thermal width (2/c^T /m) 1//2 , where T 
is the temperature at each grid point and m is the mass of the ion emitting 
the line. The spectra are integrated along a given line of sight to obtain two- 
dimensional spectral maps. The profile moments with respect to velocity are 
evaluated in order to obtain maps of line intensity, line shift, and line width. 
Based on a 30 minute time series of maps of line intensities and line shifts, the 
root-mean-square (rms) fluctuations at each spatial pixel within the respec- 
tive map are evaluated. We investigate the correlation between the average 
Doppler shift a nd relative brightness variabilities of the emission lines as pre- 



viously done by lBrkovic et al.l (120031 ). who studied time series of spectra and 



images of quiet-Sun regions observed using the SUMER spectrometer and the 
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Fig. 1. Maps in line intensity of the transition region line Civ (1548 A) formed at 
T ps 10 5 K and the coronal line Mgx (625 A) formed at T ps 10 6 K. The panels 
show the intensity images when viewing the box from the sides, i.e., along the x and 
y axis. 



Coronal Diagnostic Spectrometer (CDS) (Harisson et al. 1 19951 ) onboard the 
SO HO spacecraft. 



3 Results 



Figure 1 shows the emissivities of Civ (1548 A) and Mgx (625 A) in 
the box integrated along the x- and y-direction 30 minutes after the start 
of the simulation. This would be equivalent to an observation at the limb. 
The emissivities are plotted on a logarithmic scale according to 0.1-10 times 
the average intensity (white: emission, black: no emission). Figure 2 shows 
the emissivities of Civ (1548 A) and Mgx (625 A) integrated along the z- 
direction, i.e. when looking at the box from straight above, and the spatial 
maps in Doppler shift after 30 minutes of simulation time. This corresponds 
to an observation near disk center. The Doppler maps are scaled from -20 km 
s _1 (red; downflows) to +20 km s -1 (blue; upflows). 

Figure 3 displays the spatially averaged values of the rms intensity fluctua- 
tions (panel a) and the the spatially averaged values of the root-mean-square 
(rms) Doppler shifts (panel b) of the 30 minute time series as a function of 
line formation temperature for the extreme-ultraviolet emission lines from Ta- 
ble 1. Also shown are the time- averaged Doppler shifts as a function of line 
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Fig. 2. Maps of line intensity and Doppler shift of Civ (1548 A) and Mgx (625 A) 
when looking at the box from straight above. 




Fig. 3. rms fluctuations of intensity (a) and Doppler shift (b) as a function of line 
formation temperature for a 30 minute time-series of spectra (solid lines: simulation 
data, dashed lines: observational results). The error bars in (a) and (b) represent 
the root-mean-square deviation of the spatially averaged rms intensities and rms 
Doppler shifts, (c) shows the Doppler shifts for each line averaged over time. The 
error bars indicate the standard deviation of the time-averaged Doppler shifts. 
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formation temperature of the emission lines (panel c). 



The signatures of our model corona (solid lines) with a higher magnetic 
flux density at photospheric levels than the Sun differ significantly from the 
signatures on the qu iet Sun (dashed lines; observed data in Fig 3a and 3b from 
Brkovic et al. 120031 . observed data in Fig. 3c from Peter & Judge Il999l ). The 
quiet Sun rms intensity fluctuations show a clear maximum around 3 x 10 5 
K (Fig. 3a). The numerical experiment with a higher photospheric magnetic 
flux density than on the quiet Sun shows fluctuations that are less strong at 
higher temperatures and a maximum at lower temperatures than the observed 
solar data. The rms Doppler shift fluctuation as a function of line formation 
temperature (Fig. 3b) shows a monotonic increase for both, the simulated data 
(solid line) and the observed data (dashed line). However, the increase of the 
rms Doppler shift for the simulated data is less strong than for the observed 
data on the quiet Sun. The time-averaged Doppler shifts of the simulation 
increase with line formation temperature up to 7 km s _1 at logT « 10 6 K, 
whereas on the quiet Sun, a steady increase of line redshift up to 10 km s _1 
at logT=5.2 is observed (Fig. 3c). We take this as an indication that stars 
with higher activity levels than quiet Sun regions show a maximum redshift 
at higher temperatures as compared to the quiet Sun. However, lacking an 
observed Sun-as-a-star spectrum one has to act with caution when com parin g 
absolute line shifts of the Sun and solar-like stars. According to Peter (120061 ). 
the observed disk-center redshifts for the Sun exhibit a net line shift in the 
middle transition region three times higher as compared to full-disk values. 
Peter's (2006) study also revealed that the net line shift of solar transition 
regions lines is at least a factor of four lower than the net shifts found on the 
solar-like star a Cen A, probably indicating a higher magnetic activity in the 
chromospheric network of a Cen A. 

The trend of redshift with line formation temperature of lat e-typ e (dwarf) 
stars has b een investigate d in many works (e .g. W ood et al. 119971 . Redfield 
et al. 12001 Sim & Jordan |2002j) . Wood et al. (ll997T ) found that emission line 
redshifts of a Cen A and B increase with line formation temperature up to 
logT=5.2 almost indentical to those observed for the Sun and Procyon. Above 
logT=5.2, the redshifts decrease for the Sun, a Cen A, and a Cen B. However, 
for Procyon, they continue to increase. Redfield et al. (120021 ) investigated O VI 
lines formed at logT=5.5 for seven late-type dwarf stars including Procyon, 
AB Dor, a Cen A, a Cen B, e Eri, and AU Mic. Most stars show an increase 
in velocity offset with increasing line formation temperature and sooner or 
later a decrease with increasing t empe rature, the same general behavior as for 
solar line shifts (Peter & Judge. Il999l ). However, the O VI line shift in a Cen 
A and e Eri is highly redshifted with respect to other high-temperature lines, 
such as O V. For the cooler stars in their sample the O VI redshifts are found 
to be close to zero. At lea st for a Cen A this is a res ult in contradiction to the 
one of Wood et al. (119971 ). In Redfield et al. (120021 ) transition region lines of 
stars like a Cen A and Procyon do not suggest a drop-off in redshift over the 
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observed temperature range (up to logT=5.5; Redfield et al. 120021 ) as found for 
solar data (maximum redshift around logT=5.2), a trend consistent with the 
model results shown in Fig 3. As we do not explore a range of magnetic fluxes 
at the stellar surface, we cannot comment on the connection of e.g. magnetic 
activity and temperature or emission measure in stellar coronae yet. 

Table 1 

Wavelengths and line formation temperatures of emission lines synthesized in this 
study . 



line 


wavelength |AJ 


log(T/[K]) 


Sill 


1533.4 


4.44 


Si iv 


1393.8 


4.88 


Cn 


1334.5 


4.59 


Cm 


977.0 


4.91 


Civ 


1548.2 


5.02 


Oiv 


1401.2 


5.18 


Ov 


629.7 


5.38 


Ovi 


1031.9 


5.48 


Ne viii 


770.4 


5.80 


Mgx 


624.9 


6.05 


Mgx 


609.8 


6.05 


Fe viii 


185.1 


5.63 


Fex 


184.5 


6.99 


Fexi 


188.2 


6.07 


Fexn 


195.1 


6.14 


Fexv 


284.2 


6.32 



4 Conclusions 



We present results of extreme-ultraviolet emission line spectra synthesized 
from three-dimensional magnetohydrodynamic models of the corona with a 
photospheric magnetic flux density consisting of a quiet Sun part with a mag- 
netic flux density four times larger than on the quiet Sun and an active region 
part which involves a spatially downscaled magnetogram. We show average 
properties of the synthetic spectra and compare them to observed properties 
of solar line profiles. 
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By investigating maps of intensity and Doppler shift when integrating along 
the line of sight through the box, one can study the morphology of the transi- 
tion region and corona. Here, we show images of the transition region line C IV 
(1548 A) formed at T m 10 5 K and the coronal line Mgx (625 A) formed at 
T 10 6 K. The coronal emission appears diffuse and continuous in the upper 
part of the box, whereas one finds cool low-lying loops concentrated in the 
transition region. This result matches observation of the quiet Sun, where the 
transition region emission is confined to the chromospheric magnetic network, 
while the emission from the hot corona appears diffuse. The dopplergram of 
the transition region line is found to be highly-structured with a much smaller 
amplitude in redshift than in the corona. These findings agree well with so- 
lar observations. They support the idea that the coronal heating process is 
due to the braiding of magnetic field lines and results from photospheric foot- 
point motions. The Doppler shifts follow from the driving of the corona by 
the footpoint motion of the magnetic field lines that has been implemented 
in the magnetohydrodynamic model. As an effect of increased magnetic ac- 
tivity in the numerical simulatio n, the maximum redshift is shifted towards 
higher temperatures. Peter et al. (120041 ) showed earlier that three-dimensional 
numerical simulations can reproduce the observed variation of line shifts with 
temperature on the Sun. Our numerical experiment with increased magnetic 
complexity indicates that stars with higher levels of activity might show max- 
imum redshifts at higher temperatures. Unfortunately, observational results 
on redshifts of stars at various levels of activity are not conclusive yet. 
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